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It is shown that for superconductors with circumference 2nR approaching BCS coherence length £ minimal 
period of the response of all thermodynamic quantities to external magnetic field is set by hc/e i.e. twice the 
corresponding value for the bulk case. This is explained by the dependence of internal energy of Cooper pairs 
on their center of mass motion which leads, in particular, to a transition offset between different current-carrying 
states. Explicit calculation of the transition offset is done for the case of s-wave superconducting cylinder with 
R » £o and turns out to be exponentially small. A possible enhancement of the effect for nodal superconductors 
is suggested. Similar conclusions should also apply to the response of charged or neutral superfluids to rotation. 



It is well known that response of superconductors to elec- 
tromagnetic field and inertial perturbations such as rotation in- 
volves such characteristics of Cooper pairs as their mass 2m e 
and charge 2e with m e and e being essentially, apart from tiny 
relativistic corrections, electron's bare mass and charge. Al- 
though for all known superconductors the interparticle spac- 
ing is actually smaller (usually much smaller) than the size of 
the pairs, notions of pair's mass and charge have good heuris- 
tic value in constructing phenomenological expressions for 
various thermodynamic and transport quantities on the basis 
of their counterparts in Bose superfluids by a simple replace- 
ment m — > 2m and e — > 2e. 

The implications of this prescription on a quantum of su- 
perfluid circulation have been experimentally confirmed on a 
variety of physical systems where internal degrees of freedom 
of Cooper pairs are either absent or irrelevant. As an example 
of the electromagnetic response one can quote classical flux 
quantization experiments in both conventional flU EJ and d- 
wave superconductors [3 1; inertial properties of Cooper pairs 
were tested, in particular, in observation of quantized circula- 
tion in 3 He-B |U[5) as well as in high precision measurements 
of magnetic field induced in a superconductor by rotation |6). 
While there are no known charged Bose superfluids so that the 
unit of flux quantization established in [ 1 , 2 , 3 1, namely hc/2e, 
does not have its Bose counterpart, experiments with rotating 
superconductors [6| do have their analog-the Hess-Fairbank 
effect-in rotating 4 He and seem to well establish the circula- 
tion unit of A/mHe for 4 He and h/2m e for niobium in 
Ref.|6]thus confirming to the prescription indicated above. 

On the theoretical side flux quantization in superconduct- 
ing samples of annular geometry has been attributed to the 
condensation of pairs of particles into states with different 
momenta of the center of mass of the pairs J9j QU). In ther- 
modynamic equilibrium, as the external field changes, dis- 
crete sequential transitions between these states lead to flux- 
periodic dependencies of all quantities characterizing the an- 
nulus, e.g. to a periodic dependence of induced magnetization. 
Even without detailed knowledge of the ground state and re- 
lying only on the gauge invariance principle one can conclude 
that the fundamental period of such dependencies will be hc/e 
||9] [TT|. However if the ground state possesses pair correla- 



* vakaryuk@uiuc.edu 



tions of the type mentioned above with all or most pairs be- 
ing in the same center of mass state then the response of the 
system will contain substantial hc/2e harmonic rendering the 
minimal flux period of the dependencies to the same value. 

It has been noted a long time ago that there is no fundamen- 
tal reason behind minimal flux periodicity being hc/2e l9l [TD . 
Indeed calculation of the response of an s-wave ring under the 
assumption of independence of relative motion of electrons 
on the flux |12| as well as calculations of Little-Parks effect 
in i-wave rings [13, 14 1 and magnetization of mesoscopic d- 
wave loops lfT31 [161 have observed absence of hc/2e periods 
in the corresponding quantities. The simple argument pre- 
sented here shows that this situation, namely doubling of the 
minimal period in the response of all thermodynamic quanti- 
ties, is quite generic and happens irrespective of the particular 
form of pairing, interactions, temperature effects etc. when the 
circumference of the superconductor becomes small. For an s- 
wave superconductor where pair wave function (WF) decays 
exponentially, the effect becomes noticeable when the circum- 
ference approaches BCS coherence length However for a 
nodal superconductor, such as af-wave, there are directions in 
the real space where decay of the pair WF is algebraic which 
suggests that the effect may be noticeable at even larger than 
£o values of the circumference. Similar analysis goes through 
for rotating superconductors where the response is periodic in 
the rotation velocity with the minimal 'unbroken' period equal 
to h/2m e . 

Basic geometry considered in this work is that of a hol- 
low cylinder with the wall thickness d being smaller than the 
London penetration depth A^ and the radius of the cylinder 
R » d; same conclusions can be reached for an arbitrary ge- 
ometry requiring only axial symmetry and d «: Al condition. 
Single particle states are specified by a set of three quantum 
numbers (m, n) with hm being a projection of angular momen- 
tum along the symmetry axis and n = {n.\,ni) describing two 
other degrees of freedom responsible for the motion along and 
perpendicular to the symmetry axis. 

In thermodynamic equilibrium periodicity in the response 
of TV paired fermions is attributed to the transitions between 
ground states of the following type 
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which describe condensation of N/2 spin-singlet pairs with 
the pair's angular momentum along the symmetry axis equal 
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FIG. 1 : Schematic behavior of the internal part of pair WF as a func- 
tion of the relative angular coordinate for even (red) and odd (blue) 
values of the pair's angular momentum. High frequency oscillations 
are not shown, (a) The circumference of the cylinder is much larger 
than the decay length of the pair WF. (b) The circumference of the 
ring is comparable with the decay length. 

to tmiQ. Pairing with respect to the other two quantum num- 
bers n is chosen in a standard way and connects state n 
with its time reversal n. Spin singlet symmetry of two- 
particle state requires that variational parameter ^ ,(mo) satisfies 
X-m+m ( n ' ™) - /tm° n )- The same condition is obeyed by 
the pair wave function in the momentum representation F (mo) 
defined as 

F^ o) (n, n) = <flw ; - m + mo l a n, m T) ( 2 ) 

where the quantum mechanical average is taken between the 
states ([TJ| with N/2 - 1 and N /2 pairs. Structure of the pair 
WF is more apparent in the coordinate representation where 
it factorizes into two parts which describe internal motion and 
the motion of the center of mass of the pair: 

F(x 1 e 1 ;x 2 e 2 ) = e^V^f^ixux^ey-Oz) (3) 

Here is a polar coordinate of 2-th particle and x t repre- 
sents two other coordinates conjugate to the quantum num- 
ber n. Pair WF must be single valued in both 6 t and 82 and 
symmetric under particles permutation X\6\ ^± xiOj- The 
former condition implies that internal part of the pair WF, 
f < - ma \x\,X2,6), must be 2n periodic (antiperiodic) for even 
(odd) parity values of mo as a function of the relative coordi- 
nate 6 = 61 - 2 - 

Similar to center of mass decomposition can also be 
written for the coordinate representation of the variational pa- 
rameter x (m) which suggests that energy expectation value 
taken on state ([1} can be factorized as well. As will be shown 
below this is indeed the case. Assuming a uniform mag- 
netic flux <p through the annulus, the general form of a spin- 
independend single particle spectrum in axially symmetric ge- 
ometry with thin walls is 

s^(m, n) = [an) + (m + Qlfof] (4) 

where 0o is the single particle flux quantum <pQ = hc/\e\ and 
£(n) is some dimensionless function. For such single-particle 



dispersion, independent of a particular form of interparticle 
interactions, the expectation value of the ground state energy 
evaluated on the state <[XJ can be represented as 

£ (mo) W = [e {mo) + (mo/2 + cp/fo) 2 ] (5) 

In the expression above dimensionless parameter e ( '" o) repre- 
sents flux-independent contribution coming from the internal 
motion of the particles in the pair. Gauge invariance demands 
same parity values of s ( " l0> to be the same which allows one 
to introduce two in general independent constants, e (0) = e (2m) 
and e (I) = s ( 2m+l \ describing even and odd parity states. 

Another contribution to the system's energy comes from the 
motion of the center of mass of the pair and is represented by 
the second term on the r.h.s. of eqn. Q. It is usually assumed 
that the two contributions are independent, namely, that the 
energy associated with the internal structure of the Cooper 
pair is independent of the motion of the center of mass of the 
pair. It is the violation of this assumption which leads to the 
breaking of hc/2e periodicity in small samples. 

To see why for small rings the internal energy of the Cooper 
pair depends on its center of mass state it is instructive to 
turn to the real space representation of the pair WF (|3]l whose 
internal part is schematically shown on fig. [TJfor X\ = X2- 
For the circumferences much larger than the characteristic de- 
cay length of the pair WF fig. []J applies and the internal en- 
ergy being independent on the overall phase of WF is practi- 
cally the same for odd and even parity states. Then the cir- 
cumference becomes comparable with the characteristic de- 
cay length, the continuity condition for the pair WF requires 
it to adjust as schematically shown on fig. [TJ? [20]. The inter- 
nal energies evaluated on such states are different. 

In the simplest case of s-wave superconductor the decay of 
the pair WF in real space is isotropic and exponential with 
the decay length equal to the BCS coherence length £o- The 
internal energy difference for this case will be evaluated below 
under assumption R » and will turn out to be exponentially 
small. However for a superconductor with a nodal structure in 
the gap there are directions in the real space where pair WF 
decays algebraically which may lead to an enhancement of 
the effect. 

Coming back to eqn. <[5j and introducing superfluid velocity 
v s and superfluid density p s for the pairs of particles which 
occupy volume D. by usual relations v s = 5^1?( m o + 
p s = m e N/Q. it becomes 

O- 1 £«(0) = ^P^ o ' 1) + ip^ (6) 

Being written in this form the expression for energy can, in 
fact, be generalized to include effects of non-zero tempera- 
ture or other pair breaking perturbations by introducing corre- 
sponding changes in the superfluid density p s as well as other 
terms describing e.g. behavior of the excited component of 
the system. Such generalizations will not change qualitatively 
conclusions reached below and the following discussion will 
be limited to the simplest case specified by eqn. ([SJ. 
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FIG. 2: Dependence of the induced flux on the external flux for rings 
of small and large sizes. 

The discussion so far has led to the expressions |5]) and (|6| 
for the ground state energy as a function of pair's angular mo- 
mentum hniQ and total flux (p. To find thermodynamically sta- 
ble value of mo which is realized at a fixed external flux <p e 
one needs to construct corresponding thermodynamic poten- 
tial which is in this case Gibbs potential G. Although general 
expression for G can be easily written down, for the purpose 
of this discussion it is enough to notice that, in thermodynamic 
equilibrium, the transition between mo - and mo = 1 states 
occurs at the following value of the external flux: 
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where e <01) is the dimensionless internal energy of the Cooper 
pair for the states mo = 0,1 defined by eqn. |5]). Appearance 
of parameter a = I + 4-Rd/A 2 ^ in (T7fc has purely geometri- 
cal reasons related to fact that, for the geometry considered 
here, field and condensate energies scale differently with R. 
As one can see from the result above, transition between dif- 
ferent center of mass states is offset relative to the bulk case 
by an amount proportional to the internal energy difference. 

Complete dependence of the induced flux obtained by min- 
imization of the Gibbs potential corresponding to the energy 
Q is illustrated on fig. [2] The effect of the non-vanishing in- 
ternal energy difference is to shift the transitions between dif- 
ferent current states by an amount proportional to +(e (I) -e (0) ) 
thus breaking hc/2e periodicity up to hc/e Oil . In an alterna- 
tive interpretation one may refer to this effect as a change in 
Cooper pair's charge. A very similar reasoning shows that the 
same change occurs for the Cooper pair mass in the rotating 
superconductors thus leaving gyromagnetic ratios like 2m e c/e 
unchanged. 

The preceding qualitative conclusions were reached with- 
out particular reference to a microscopic form of interparti- 
cle interactions. However for a quantitative analysis, which is 
most easily done in the momentum representation, a choice of 
specific microscopic model is needed. Using single-particle 
dispersion relation Q and a contact interparticle interaction 
potential V energy expectation value on a many-body ground 
state ([]} is 

E (mo) (<P) =J]2 £ ^m, n)\v ( ^ no/2 (n, nf 

-Z y ^/ 2 (-'-)Ct/ 2 ("'-) (8) 
where quantities v and T are defined through earlier intro- 
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X^'{n,n)F^ 0) \n,n) and T^ o/2 (n,n) = F% a \n,n). De- 
fined in such way v and T are symmetric with respect to m, 
e.g. v m (n, n) = u_,„(n, n) etc. In what follows, to avoid com- 
plicated notation, n dependence of all quantities will not be 
indicated explicitly unless mentioned otherwise. 

Noticing that normalization condition for state ([T]l requires 

2 Z \ v m'-l, /2\ 2 =N, the flux dependence in (8j) can be sepa- 
rated out thus establishing a connection witn the expression 
for energy <[5j used in the previous analysis. Identifying flux 
dependent part with the the center of mass contribution and 
everything else with the internal energy contribution one thus 
obtains for the latter 
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and the interaction constant is now written in units of 
h 2 /2m e R 2 . It is also apparent that e (mo) defined by the equation 
above is the same for the same parity states, as was indicated 
above. 

The actual forms of y (m °^ and <7 r(m >>'> can be found by en- 
ergy minimization in the same way as it is done in the stan- 
dard BCS treatment which leads to \v^° } \ 2 = \v m (A {ma ) n (m) )\ 2 , 
Tm a) = T m (A {m) ,^ mo) ) where, for a given pair state m , the 
gap parameter, A (m °', and the chemical potential, should 
be determined self-consistently from gap and normalization 
equations and are, in general, different for odd and even parity 
states. For fixed A (mo) and fi {mo) functions v {mo) and T im) are 
given by standard BCS expressions. 

In the momentum representation the reason for a non- 
vanishing value of the internal energy difference between odd 
and even parity states and, as a result, doubling of the flux 
periodicity, is less transparent than in the coordinate represen- 
tation but from a mathematical point of view it can be traced 
down to two main factors. In the first place, the actual val- 
ues of m arguments which v and T are being summed over 
in (|9]l are shifted by 1/2 for odd and even parity cases; sec- 
ondly, there is a difference between odd and even values of 
the gaps and chemical potentials. It can be shown that the 
latter difference being itself caused by the former in the gap 
and normalization equations can be ignored in the leading ap- 
proximation. The actual calculation of the energy difference 
can be found in the Appendix. Performing summation over m 
with the help of Poisson summation formula and integrating 
afterwards with respect to two other quantum numbers n it is, 
in fact, possible to get an analytic expression for e (1) - e (0) in 
the limit R » £q. With the relative accuracy %o/R one has 
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where kf is the Fermi wave vector and £o is defined through 
the gap parameter A as £o = fivp/nA. It should be emphasized 
that the accuracy of eqn. (TTOJ does not allow one to distin- 
guish between A (1) and A (0 Mor and // 0) ) since the differ- 
ence A (1) - A (0) is itself exponentially small. There might be 
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however a substantial difference between either of A* ' 1 ' and 
the value of the gap for the corresponding bulk material. 



One should observe that the sign of eqn. ( 10 1 is a rapidly os- 
cillating function of the radius which is due to the oscillations 
of the pair WF in real space. Similar oscillatory effects have 
been found theoretically in Ref. [ 17] for the dependence of the 
gap on the thickness of a thin superconducting film. Strong 
dependence of the sign of the transition offset on the ring's 
size would make it more difficult to observe it on an ensemble 
of rings, such as the one used in [ 18 1 for Little-Parks measure- 
ment, because of possible variations in the rings' sizes. 

Another limitation for experimental observation is the mag- 
nitude of the effect 



10 1 one may expect 
3ft ; larger 



According to eqn. 
transition offset to be of the order of 0. 1% for R 
corrections for the same value of /?/ft will most likely be 
achieved for nodal superconductors. Although reducing the 
radius will make the effect more pronounced, a little more 
care should be taken in extrapolating eqn. ( 10 1 to the ring 
sizes equal or lesser of the coherence length. This is due to 
the fact that when the diameter of the ring becomes equal to 
the coherence length the superfiuid velocity needed to screen 
a quantum of flux reaches its critical value thus leading to the 
destruction of superconductivity around cp e = <po/2 {19j . 

It is interesting to notice that, apart from the small sample 
size, there might be other ways to introduce correlations be- 
tween the center of mass and internal energies of the Cooper 
pair. One of them could be a strong spin-orbit interaction 
where the small size limitation would not apply. 

To conclude, it was suggested that in small superconduc- 
tors, due to the dependence of internal energy of Cooper pairs 
on the center of mass state, the minimal flux periodicity is 
hc/e, twice the usually attributed value. The doubling of the 
periodicity is due to the offset of the transition between differ- 
ent current states. The magnitude of such offset was calculated 
for s-wave pairing. It was also suggested that other things be- 
ing equal, the effect will be more pronounced for nodal super- 
conductors. 

The author wishes to thank Tony Leggett for guidance 
and David Ferguson and Alexey Bezryadin for useful discus- 
sions. Support from NSF grant DMR03-50842 is gratefully 
acknowledged. 



APPENDIX 

In this appendix expression ( [T0| for the internal energy dif- 
ference between odd and even parity states of a large, R » ft, 
i-wave superconducting cylinder or ring will be obtained. Un- 
less explicitly otherwise stated all quantities with dimensions 
of energy are written in units of h 2 /2m e R 2 . 

As it has been already remarked the minimization of the 
expectation value of energy given by eqn. (IBl leads to the the 
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the standard BCS expressions: 
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with e,„ = £ + m 2 - /j. being a single-particle energy counted 
from the chemical potential which corresponds to a given 
value of m.Q. 

Now, using eqn. (|9j, the internal energies for even and odd 
parity states are given by 

Nsf = 2 + '« 2 )k,,ol 2 - V J] r m , r^ (A.5) 
for even parity states, and 
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for odd parity states. The subscript or 1 in all quantities 
indicates dependence on A (0 \ f/® or A (I) , respectively. 
The internal energy difference can then be written as 
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energy (mo = state) as A and fi are varied from their mo = 
to mo = 1 values. Since el 0) is the equilibrium i.e. minimal 



value of e (0) , the last term on the r.h.s. of eqn. (A.7 1 is of the 
second order in A (1) - A (0) and, as it will be seen below, can 
be ignored relative to the first two terms in the limit R » ft. 
Thus to evaluate the internal energy difference in the leading 
order the energies themselves can be taken at the same values 
of A and At this accuracy any difference between A <0) and 
A (1 ' etc. is neglected; however, there might be a substantial 
difference between these values and the corresponding bulk 
parameters. 

Using the gap equation A = V^Tm-i/2, in the leading ap- 
proximation eqn. (A.7 1 becomes 



Nis^- e (0) ) = 2^ + (m-l/2) 2 )t, 2 _ 1/2 (A.8) 

- + m2 >i - 2A Yi {Tm - x i 2 ~ Tm) 



where v m and T m given by eqn.'s ( | A.3| > and \KA\ . The sums 
above run over all integers m as well as over two other quan- 
tum numbers n which represent other-than-azimuthal part of 
the dispersion ((n). With the help of the Poisson summa- 
tion formula the sum over m is converted to an integral; sum- 
mation over n is replaced by integration using the standard 
rule Yjn ~ * J d£,giiO where gi(£) is the (dimensionless) den- 
sity of states for f (n) and can be considered to be a constant. 
Keeping only leading exponential term in the Poisson summa- 
tion series one obtains: 



where, for given values of A and ft, v m and T m are specified by 
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where the following definitions are made 



J d£gi J dx 



e(x, £) e 



2nix 



and 



= 4 fagij ix ^e 2 (x,0 + A 2 - e(x, 



(A. 10) 



(A.ll) 



with e(x, = ( + x 2 - fi. It is intuitively plausible and can, 
in fact, be shown that ~ £o/R <s 1 so that the problem 
reduces to the calculation of the following integral: 
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where and from now on f will be counted relative to the chem- 
ical potential. The rest of the discussion is devoted to the cal- 
culation of the above integral under assumption R » £o- 

The integration over x in eqn. ( |A.12| i is transformed by clos- 
ing the integration contour in the upper half of the complex 
plane where the denominator + x 2 ) 2 + A 2 has two branch 
points. Choosing the branch cuts to point outwards away from 
the origin the integral over the real axis reduces to four inte- 
grals along both sides of the two branch cuts. The two in- 
tegrals along one branch cut are complex conjugate of those 
along the other branch cut. Introducing new integration vari- 
able t along the first quarter branch cut by z = Zo(t +1), where 
Zo is the first quarter branch point, the oscillating exponent in 



(A. 12 1 acquires a decaying part: 
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where % denotes the real part of the corresponding expression 
and the decaying exponent is given by the imaginary part of 
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The essential for the integration over t region in (A.13 1 is 
determined by the decaying exponent e~ 2mj0 ' and extends from 
t m m = to f max ~ (2m/o) _1 . As a function of f, yo is mono- 
tonically increasing in the region (-//, oo) reaching its min- 
imum at ( — -fi with the value yo(-/i) ~ A/2-07; at the 
same time xo is monotonically decreasing with the maximum 
value xo(-fj) ss yfji ll22l . Restoring dimensions i.e. supplying 
h 2 /2m e R 2 denominators to A and fi gives 



yo(-fi) ~ R/ni-o, x (-fi) ~ k F R 



(A.16) 



where £o is the BCS coherence length £o = nv F /nA and u 2 = 
2jd/m e [23]. 



The above analysis suggests that, under condition R » £o, 
the integrand in ( |A.13[ ) can be significantly simplified because 
in the important for the integration region f « 1 holds. After 
some algebraic manipulations followed by a Wick rotation of 
t variable to make the algebraic part of the integrand real for 
£ < 0, the expression for the integral takes the following form 
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with further corrections being suppressed by a factor of order 
of ft//?. 

The integration over £ is considered separately for intervals 
of positive and negative values of After rescaling variable 
t by A/|£| in each interval the non-trivial dependence of the 
integrand on A and ( is gathered in the exponent. For the 
integral over f > region the exponent has a sharp maximum 
and the integral can be evaluated by the saddle point method 
leading to its magnitude being of order of exp {—2n VA). As 
it will be seen below the integral over negative values of f 
is of order of exp (— 2ttA/ yfjl). Since the ratio of the the two 
exponents is y/A/fi <K 1 the integration over f in ( A.17 1 can 
be limited to the interval ( < leading to 
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where 



(A.19) 



X(i,t) = 2ny (O + MO^t 
Y({,t) = -2kx (O + yoiO^t (A.20) 



with xq and yo given by eqn.'s ( A. 14 1 and ( A.15 1. 

Expressions for X((, t) and Y{(, t) can be significantly sim- 
plified in the limit \(\ » A. In particular for f = -yu 



where R is defined as 
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and in dimensional units is equal to R/^o, the ratio of the ra- 
dius of the annulus to the coherence length. As it has been 
already mentioned it is assumed that R » 1. 

The function in the exponent in eqn. ( |A.18| >, X((, t), is a 
positive monotonically growing function of f for f < which 
reaches its minimum at the lower integration limit ( — —jj, and 
diverges at the upper limit f = 0. To take the advantage of 
the simplified forms of X and Y at f = -jj, one can perform 
repeated integration by parts in the integral over f obtaining 
in that way following asymptotic expansion: 
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where the prime sign denotes differentiation with respect to 
( and function Z((, f) is defined as the £ dependent non- 
exponential part of the integrand in ( |A. 1 8| >: 



Z(CJ)= — sinF(<r,/) 



(A.25) 



Despite complicated at first sight t dependence of the inte- 
grand in ( |A.24[ ) the integration over t can be carried out notic- 
ing that, because of the decaying exponent, the interval rel- 
evant to the integration is limited by R l <K 1. Dropping t 
dependence relative to constants of order 1 and integrating the 
expansion term by term gives with R l relative accuracy 
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where S (jj.) is defined as 
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X>({,0) 



(A.26) 



(A.27) 



The sum above cannot be limited to a finite number of terms 
because the n-th term in the sum is of order of [T ' 2 (ji/A) . 
However by rearranging summation one can notice that S (ju) 
satisfies following differential equation 



S(ji)=Z(-fi,0)- 



d S(jj) 
dp X'(-n, 0) 



(A.28) 



or, given explicit form of X(-/j, 0) and Z(-/z, 0), eqn.'s (A. 21 



and ( A.25 1, the differential equation for S (ji) becomes 



S Qx) = sin In ^ - ^ [// 3/2 S Qx)] 



(A.29) 



It can be checked by a direct substitution that up to the terms 
of order of Aj/i this equation has the following solution: 



.3/2 



! 2n yfjl 



(A.30) 



so that the expression for the integral I reduces to 

/ = ^3/2^1/2 e -2R cos27T ^l (AJ1) 



where g 2 is the density of states at the chemical potential 
for other-than-azimuthal part of the dispersion (,{ri). Using 
single-particle dispersion law for a thin walled cylinder and 
reverting to dimensional units one can see that g 2 /i 3 ^ 2 = \N 
so that 



3./V 



1/2 -2R/h 



cos InkfR 



(A.32) 



Comparing this equation with eqn. ( |A.9| l gives the result ( |T0] > 
for the internal energy difference quoted in the main text. 
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